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Masses of the r] and r\' mesons are estimated in Nf ■ = 2 + 1 lattice QCD with the non- 
perturbatively 0(a) improved Wilson quark action and the Iwasaki RG-improved gluon action, 
using CP-PACS/JLQCD configurations on a 16 3 x 32 lattice at j3 = 1.83 (lattice spacing is 0.122 
fm). We apply a stochastic noise estimator technique combined with smearing method to evaluate 
correlators among flavor SU (2) singlet pseudoscalar operators and strange pseudoscalar operators 
for 10 combinations of up/down and strange quark masses. The correlator matrix is then diago- 
nalized to identify signals for mass eigenstates. Masses of the ground state and the first excited 
state extrapolated to the physical point are m^ = 0.545(16) GeV and m„i = 0.871(46) GeV, being 
close to the experimental values of the rj and rj' masses. 
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1. Introduction 

The U(l) problem [jl|] is an outstanding issue in hadron spectroscopy. The large mass of the 
tj' relative to the % is believed to be related to QCD vacuum structure and anomaly of the axial 
current. In order to reproduce the tj' mass from first principles calculations, a number of lattice 
QCD simulations have been carried out. Early attempts were made in quenched QCD. Recent 
studies [Q, ||, ^] shifted to full QCD with two degenerate up and down (ujd) quarks. 

In the real world the role of strange (s) quark is important since the mass difference between 
the s quark and ujd quarks is the cause of mixing between the singlet and octet pseudoscalars, 
that leads to the physical tj and tj' mesons. In this article, we report on an attempt to estimate 
masses of the tj and tj' in Nf ■ = 2 + 1 QCD, in which degenerate ujd quarks and an s quark are 
treated dynamically. The mixing effect is taken into account by diagonalizing a correlator matrix 
for operators coupled to both tj and tj'. We use Wilson quark formulation. An attempt to calculate 
the tj' mass with staggered fermions is made in Ref. [0]. 

The calculation is carried out on a set of gauge configurations generated previously for a study 
of the flavor non-singlet hadron spectrum and light quark masses in Nj ■ = 2 + 1 QCD [JSj] . This report 
presents results for the tj and tj' masses calculated for the coarsest of the three lattice spacings. 

A technical difficulty in the study of the tj and tj' masses lies in the evaluation of the two-loop 
quark diagram contribution to the singlet part of propagators. Due to configuration fluctuations in 
the two-loop diagram, errors of the tj and tj' propagators increase rapidly as time slice. On the 
other hand, we have to suppress contamination from higher excited states. Calculations have to 
be setup in a way that signals of the tj and tj' appear at small time slices. For this purpose, we 
combine a smearing method with a stochastic noise estimator technique (SET), as carried out in 
Ref. [§. 

The organization of this article is as follows. Numerical calculations including our procedure 
to calculate the two-loop diagram are described in Sec. ||. In Sec. |5| we give an account of analysis 
procedure and the tj and tj' masses at the physical point. Conclusions are given in sec. ^ Data 
presented here should be regarded as preliminary since detailed analysis is still ongoing. 

2. Numerical Calculations 
2.1 Configuration 

We use a set of Nf = 2 + 1 full QCD configurations ^ generated with the Iwasaki RG- 
improved gauge action and the non-perturbatively 0(a) improved Wilson fermion action at j3 = 
1.83 on a 16 x 32 lattice. Polynomial hybrid Monte Carlo simulations have been carried out for 10 
quark mass combinations; we choose five hopping parameters K uc i = 0.13655, 0.13710, 0.13760, 
0.13800 and 0.13825 for the u/d quarks, and two K s = 0.13710 and 0.13760 for the s quark. The 
pseudoscalar to vector mass ratio ranges from 0.78 to 0.61. The lattice spacing is a = 0.122 fm 
which is determined from experimental values of the %, p and K meson masses. 

The trajectory length is 7000 - 8600 for each mass combination. Measurements are carried 
out at every 10 trajectories. Errors are estimated by the jack-knife method with a bin size of 10 
configurations (100 trajectories). 
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Figure 1: Error of the two-loop diagram versus number of noise sets for (K s = 0. 13710, A^/ = 0.13825). 
The vertical scale is arbitrary. 



2.2 Two-loop quark diagram 

When one determines the tj and tj' masses from propagators around a time slice, contamina- 
tion from higher excited states has to be small at that time slice. For this condition to be satisfied, 
it is desirable that contamination is small also for octet pseudoscalars. In Ref. [^Jj, we adopted an 
exponential-like smearing function 

K(n-m) =Aexp(-B|n-m|) for n^fh, K0) = l, (2.1) 

with A = 1.2 and B = 0.1, and find that effective masses of octets obtained with doubly smeared 
source and unsmeared (point) sink operators reach an approximate plateau at a small time slice; 
the mass estimated at t m 4 is only a few percent different from accurate value of the ground state 
mass. With this observation in mind, we construct the two-loop diagram using a doubly smeared 



loop diagram at the source with the kernel ( |2.1| ) and an unsmeared loop diagram at the sink, and 
combine it with the corresponding one-loop diagram. 

The two-loop diagram is estimated by a combined method of SET and smearing applied to 
configurations fixed to Coulomb gauge. We prepare a U{\) random number exp(/0(H,?)) for each 
site (n,t) of the whole lattice and solve a quark propagator q(n,t) for the random source. The 
propagator is then smeared twice at a source time slice f src , multiplied by the inverse of the random 
number and summed over sites on the time slice, 



(2.2) 



£ exp(-id(lj STC ))K(l-m)K(m-n)q(n,t SIC ). 

l,rh,H 

The loop diagram £(f sn k) at a sink ti me slice 4nk is calculated similarly without smearing. The 
two-loop diagram is estimated by averaging the product of the two loops over 20 stochastic noise 
sets generated for each color and spin at a time, 

1 



G2-loop(/) 



L 3 r 



£ (Tr(75Ctosn k ))Tr(7 5 ^(4rc)))r 



(2.3) 



In any stochastic noise methods, error consists of stochastic noise and configuration fluctua- 
tions. In our case, the error of the two-loop diagram decreases rapidly with a number of noise sets 
and reaches a plateau for a range of time slices we use for later analysis. See Fig. |I[ This means 
that our final error obtained with full 20 noise sets is dominated by configuration fluctuations. 
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Figure 2: Diagonalized propagator for (K s = Figure 3: Effective masses of the 77 and 77' for (K s = 
0.13710, K ud = 0.13825). For off-diagonal part, 0.13710,^, = 0.13825). Those for % and r/.„ are 
D l2 {t) in eq. ([O]) are plotted. overlaid. 



3. Analysis and Results 

3.1 t] and rj' masses at simulation points 

In order to derive mases of the r\ and r\' , we consider flavor SU (2) singlet and strange pseu- 
doscalar operators 

r] n = {uy 5 u + dy 5 d)/V2, r} s = {sy 5 s), (3.1) 



and calculate the 2x2 correlators 



G(t) 



<M(0) <M(0) \ 

nf(Oi» 5 (o) (0»?/(o) J ' 



(3.2) 



where superscripts S and P denote smeared and point operators. If the correlator matrix is domi- 
nated by the contribution of the ground state (rj) and the first excited state (r)') at a time slice t, we 
obtain a relation 



(CG{t)G-\t )C- 



8ij exp(-m,-(f - 1 )) i, j = rj , 77', 



(3.3) 



where the matrix C diagonalizes the basis operators ri[ (k = n,s) to give the physical ones Oi 
(i = 7] , 7]'), Oj = Y*k=n,s CikTIk* an ^ f o is a reference time slice. The matrix C is independent of time 
slice. Selecting a representative time slice tjj, we estimate the matrix C by diagonalizing the matrix 
G{t D )G~ l (to). We then calculate the "diagonalized" propagator 



D ij (t) = (CG(t)G-\t )C- l ) i 



(3.4) 



Note that Dj,-(fo) = 8u and D;,(f D )=0 (/ 7^ 7) by definitions. 

Figure ^| shows the diagonalized propagator obtained with ?o = 4 and tp = 3 for the mass 
combination of the heaviest s-quark and the lightest w/J-quarks. Diagonal parts (signals for tj and 
Tj') decay exponentially as indicated by effective masses plotted in Fig. |3[ while off-diagonal parts 
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Figure 4: Effective masses of the 77 (open symbols) and rj' (filled symbols) for K s =0.13710 (top panels) 
and 0. 13760 (bottom panels), u/d-quark mass decreases from left panels to right panels. 



are negligible for all t. In other words, the relation (3.3) is realized. This result suggests that the r\ 
and v\' contribution dominates our correlator matrix at small time slice, e.g., already at t » 2. 

In Fig. ||, we overlay effective masses of an octet pseudoscalar ujsd (71) and those of syss 
determined by ignoring the two-loop diagram contribution ("ri ss "). The mass of the v\' is much 
larger than masses of the n and 7] ss , which demonstrates the importance of the two-loop diagram 
contribution for the large mass of the 77 ' 

Figure [| shows effective masses of the Y\ and r]' for all quark mass combinations. Because 
effective masses of the rj show approximate plateaus around t rj 4, are determined from fits 
in the range t = 3 — 5. For v\' , plateaus are clear for the lightest w/<i-quarks, while they are not 
so apparent for heavier w/<i-quarks. Nonetheless off-diagonal parts of the diagonalized propagator 
are consistent with zero already at t = 2 for all cases, which suggests that contribution from higher 
excited states is small. Therefore, we estimate m^i from fits in the range t = 2 — 4. Table | 
summarizes masses of the v\ and v\' . 

For degenerate u/d and s quarks, the T] is a pure octet state, and therefore should agree 
with m^. Values of in Table [T] are determined by the procedure explained above even for these 
cases. On the other hand, one-mass fits to T] ss propagators in the range t = 8 — 16 give m^ st = 
0.6251(13) and 0.5136(09) for K ud = K s = 0. 13710 and 0.13760 respectively. They are larger than 
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Table 1: Masses of the 77 and rj' mesons. 
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Figure 5: Pseudoscalar meson masses versus 1 /K uc ( for /f v =0.13710 (left figure) and 0.13760 (right figure). 
Stars marked at the physical u/d quark mass point indicate experimental values. 



m-q by about 2 percent. We regard the 2 percent as an estimate of systematic error in arising 
from excited state contaminations. 

3.2 Chiral extrapolations 

Chiral fits to and m^i are made using low order polynomial functions of quark masses de- 
fined by m ud = (1 /K ud — 1 /K c ) /2 and m s = (l/K s — 1 /K c ) /2 (K c is the critical hopping parameter). 
Masses of pseudoscalar mesons are plotted in Fig. Q as a function of \/K U( j. The figure suggests 
that m^i shows a linear dependence in m ud within our errors, so we assume a simple linear function 

"V =A + B ud m ud + B s m s . (3.5) 

We obtain # 2 /dof = 1.4, which is acceptable. We include a quadratic term of m u d for m^, 

m n = C + D ud m ud + D s m s + Em 2 ud , (3.6) 

because a linear function does not fit to data due to curvatures shown in Fig. |5[ The quadratic term 
reduces % 2 of the fit greatly to # 2 /dof = 1.0. We note that linear fits to m 2 rjl give masses at the 
physical point consistent with those from the fits above within 3a for the rj and la for the rj'. 
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Masses of the v\ and v\' extrapolated to the physical point, as determined from experimental 
m n , m p and nix, read 

= 0.545(16) GeV, (3.7) 
«V = 0.871(46) GeV. (3.8) 

4. Conclusions and Future Plan 

We have estimated the masses of the v\ and tj' in Nf = 2 + 1 lattice QCD, albeit at one lattice 
spacing. Our result m,j=0.545(16) GeV is consistent with the experimental value of 0.550 GeV. The 
tj' meson mass m,j/=0.871(46) GeV is also close to experiment, 0.960 GeV. This is encouraging 
since the small difference of 100 MeV between our value for m^i and experiment may well be 
accounted for by systematic errors. 

Finite lattice spacing is one of possible origins of systematic errors. Another shortcoming of 
this work is a relatively large value of u/d quark masses. We expect to improve our program for 
the U(l) problem using CP-PACS/JLQCD configurations already generated for finer lattices and 
configurations which will be generated by the PACS-CS collaboration [^] soon for lighter quarks. 
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